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Summary of DE 195 30 646 CI 

DE 195 30 646 CI discloses a learning method for a recurrent 
neural network. By the application of this method, the 
returns of recurrent neural networks are splitted. This 
provides additional inputs receiving training data. The 
training data are obtained on the basis of error 
distributions and the application of a Monte Carlo method. 

DE 195 30 646 CI represents technological background with 
regard to the use of neural networks. It does not disclose a 
method for detecting the modes of a dynamic system with a 
drift segmentation model as claimed in the above U.S. patent 
application. 



PAGE BLANK Iusfto) 



® HUNDESREPUBL1K 62) Patentschrift 

DEUT3CHLAND @ Q£ 19530646 CI 




DEUTSCHES 
PATENTAMT 



Int CL»: 

G 06 F 15/18 

G05B 13/02 

o 



© Aktsnzsichen: 19530848.5-53 

@ Anmeldetag: 21, 8.95 

® Offenlegungstag: — 

© Veroffentiichungstag 

der Patentarteiiung: 17. 10. 96 



Innerhafb von 3 Monaten nach Veroffentltehung der Erteilung kann Snspruch erhoben warden 



I Patentmhaber: 

Siemens AG, 80333 Munchen, DE 



Erf? n dsn 

Tresp, VoJker, Dr., 80539 Munchen, DE 

Fur die Beurteilung der Patentfahfgkeft 
In Betracftt gezogene Druckschriftsn: 

US 51 82 794 
TRESP, Volker, etal,: Training Neural Networks with 
Deficient Data. \n: Cowan, J.D.etaL(Hr$g.) Advances 
in Neural Information Processing Systems 6, Morgan 
Kaufman, 18S4 f S. J 28-135; 
J P 06-301 663 A,. In: Patent Abstracts of Japan 
(CD-ROM); 



<§) Larnverfahren fur ein rekurrantes neuronals Netz 

© Bel Anwendung des arfindunflflgernS&en Varfahrens war- 
den die ROckfOhrungen rakurrenter neuronalar Natza aufge- 
trennt. Fur die hierdureh gawonnenen zusdtzllchan Elngange 
am Netz warden TraEningadatan fcunstfich bereitgasteift, 
indent anhand dar baksnnten Mefiwerta aua dor zu Trai- 
ning szwecken zugafQhrtan Zaftreiha und daran bekanntar 
oder vorgegabener Fahtarvarteflungadichta fur die zusatdi- 
chen Werto FehlarvorteJIungsdlchtan berechnet warden. 
Hlerzu eind gem5& der Moote-Carlo-Methode aua dteser 
FehJaivertailungsdichte Proban zu zlehan. Mrt denan wlrd (a 
ein vorherzusagender Wart der Zaltrefhe barachnet und 
deren Mfttalwert wJrd fOr einen vorherzuaagenden Wert 
ein gesatzt. Hnaatzgablete dar HrHndung Began auf alien 
bekannten Bnsatzg ebieten neuronalar Netze. 



to 

o 

« 

to 

CD 




Q 



BUNDESDRUCKEREI 08. SS 902142/335 



IS 



DE 195 30 646 CI 



Bescfcre&ung 



Die Eifindungbezieht rich auf 

Neuronale Netze fmdea in (fie vtdfikigsten tcchnisciien GeMete Eingang. Oberall don, wo es gut aus 
komplexen tectaischen ZusammenhSngen and aus mzurdchernten Infbnnationen Entacheidungen abzulciten, 
erweisen rich neuronale Netze mis besonders gedgnet Zur Bifdung cincr odor mehrerer AnsgaiigsgrdBea 
werden dem neuronalen Netz bdspielsweise erne odcr melirere EingangsgrdBen zugefOhrt Hierztt wad em 
jofchci Netr nmadurt fflr drr qf^> n *" ¥5ngat*fiiff trainiert aasAGeflend gencraBsicrt and dananh wnd es nut 
rinem snderen Dfltnmrr ah dm ^«fr g*rf«*«i gafrfiert. Neuronale Netze erweisen rich for role I3ns at rf fill e 
abbeaondersgeeigiietdaMtnihre^ _ 

Em hfefig auftretendes Problem im Znsammenhang mit dem Einsatz von neuronalen Netzen best*t afler- 
<fingg darin, daB hiufig die Eingangsdaten zom Training, oder beim Betrieb des Netzes mcht voUstaiidig smd. 
DiSer Sachvcrhait und auch die Tatsache, daB die MeBwerte fQr den Aufbau einer Zeitreaie, wdche don 
neireraaleaNetz zugefOhrt wW 

ggbmssc der Netze crzielt wcrdciL Bd reknrrentBa neuronalen Netzen tritt insbesondere das Problem aut daB 
wegen der auf Eingange rflckgekoppdten Ausgange beim TVaining nicht aile &nginge mit Traimn^datea 
beschxekt werden kdimen. Bislang gibt es keine Methoden, welche cfieser besondercn Problematik Redmong 



Aus der Vertffentlicfaiing von TYesp, Volfcer et at: "Training Neural Networks with Deficient Data * Ins 
Cowan, JJX et sl(Hrsg): Advances in Neural Information Processing Systems 6, Morgan Kaufman, 1W4> S. 
12&— 135,sind die Gnrndlagen zum Training neuronaler Netze mh Hiife unvollstSndiger Datea be kamrt . Vor 
alien Dingen werden dabei L6sungen bevorzugt. die eine gewichtete Integration fiber die unb e karmt e n oder 
unbestnnmten EingSnge erfordem. 

Aus der US-Patentschrift US 5 182 794 ist em Lernsystem fQr rekurreotc neuronale Netze bekarmt. Dieses 
Lernverf ahren betrifft neuronale Netze. welcfae vcrstecktc Eingangs- und Acsgangsneuronen haben und fadkn- 
Bert WtchtungsfeMer fiber eine beschrSnkte Anzahl von Fortpflanznngen liber das Netzwcrk- Mit diesera 
Verfahren kSnnen gangigt: Lernmengen, weldie fQr vorwartspropagierende Netzwerke vcrwendet werden* 
uidifBrrdcurrtnteNetzweto # 

Ans der japanischen VerfiffeniBchung JP 06-301 663 A. Im intent Abstracts of Japan ist era Lernsystem fur 
seriell verbundene rekurrente neuronale Netze bekannt Dabei wird der aktive Ausgangswert ernes Aiisgangs- 
neurons 4 als Lernsignal far dnen Ansgangsznstand verwendet, und in einem fortgeschrrttenen Stadium wird 
beim Updaten der Verbindong der Ausgangswert des Neurons 4 und der Fingmi g s wertdes Lemsagnals mtt einer 

vwywhtvhnt^wwwn f^wirii trmffrfrmfaion gewkhtet . , 

Die der Erfindung zugruudeliegende Anfgabe besteht deshalb darin, em Lernverfahren anzugeben, nnt dem 
der Lenxvorgang bcimTrainiug ernes reknrrenten neuronalen Netzes verbessert werden kamL 
Diese Aurgabe wird gema^denMerionalendes Patentanytprc^ 1 gel^t 

Wdterbadungcn der Erfindung ergeben sidi ans den abhSngigea Ans [irQehen. ^ 

Bin besonderer Vortdl des erfindungsgem&Ben Verfabrens besteht dann, daB hier erstmals em Verfaiffen 
vorgestdh wird, not welchem mehr Einginge fQr das Training rekurrenter neuronaler Netzweite znr VerfB^ 
gung stehen. Durch das Verfahren werden namlich die Rfickffihrungen der rekurrCTten Netze aafeetrquit und 
dieseairfgetrenntenVeibrndnn 

se steHt das erfindungsgcml8e Verfahren hochwerttge kfinstfidi erzeugte Werte fflr cfie Zeitrdhe tsxjt Verffl- 
gung, wddie mit Hilfe der Fehlerverteihmgsdtchte und der restGdien bekaimten Werte der Zeitrdhe mit 
wekhem das Netz trainiert werdox soli beredmet werdra. , . 

Besonders vortdlhaftkann das erfindungsgerrtBe Yerfafaren such angewendetwerden.wennmehrere ndsen- 
einanderfiegende Werte dCTZeitT^^ 

Verbindungen aufeetramt werden. Zar Beredmnng dteser Werte wird vorzugsweise em fteralionsveifahren 
ang^ebeffl, bd wdchem ahnfich vorgegangen wird, wie bdder Berechnung eines Euaehrates, jedoch wird 
dann dieser Wert als beksumt vorausgesetzt und der weitere unbetonnte Wert mit diesem beredmet Die 
Iteration kann solange dnrchgefQhrt werden, bis eine hinreicheiide Genauigkeit dex zu enmttdnden Werte 
gdundenwuide, oder to das Netz mit dragewunschtenLera ^ 

Ein besonderer Vorteil des erfindtmgsgemlBen Verfahrens besteht darin, daB ausgenutzt wird, daB die 
fehlenden Werte oder die verrauschten Werte, weiehe den neuronalen Netz zugefflhrt werden soUen. Bestand- 
tdl dner Abfolge von Werten in der Zeitrdhe sind VorteOhaft kann die bekaimte FehlervertdfangSTOhrsdieui- 
lidikdt der restlichen Werte dazu benutzt werden, urn ttr den fehlenden Wert nach dem erfindungsgemSten 
Verfahren eine erwartete Fdilervertdlung imd darmt den erwartetea Wert berechnenzu konnm. 

GQnstigerweise konnen nadi dem erfindungsgemaBcn Verfahren auch fehlende Werte, die in der Zeitrdhe 
benachbart sind, bestimmt werden. DafQr ist dn iterativer Vorgang vorgesehen, der emmal den dnen Wert 
berechnet und danach den anderen mh den aus dem einen Wert gewonnenen Daten ermhtdt Gfinsti^rweise 
kann cfieser Itcrafjonsvorgang audi mehrfach durcfagefOhrt werden, damit eine hinreidiaide Genauigkeit der zu 
bestimmenden Werte gewahrleirtet isL m m .it 

Besonders kann nach dem erfindungsgemSBen Verfahren auch dn neuronales Netz trainiert werden, wdenes 
die Zeitrdhe nachbildett soil, da dabei vortdlhaft die Lernschrittweite auf die Zahl der gezogenen Monte Carlo 
Proben bezogen wird . 9 „ . . . . n t „ 

GOostigerweise wird fur das Fehlerverhalten der Emgangsgr6fien der Zeitreihe erne GauBverteihmg ange- 
nommen oder festgdegt da cfies eine Verteilung ta, welche praxisnahen Werten wdtestgehend entspricfat 

Vorteflhaft werden dnfache numerische Verfahren beim erfindungsgemSSen Verfahren eingesetzt, um feh- 
lende Werte bestimmen zu konnen bzw. um einen kflnftigen Wert in der Zeitrdhe mit fehlenden Werten 



r 



DE 195 30 646 CI 

vorhersagen zu kdnnen. 

Vorteilhaft werden nach dem erfindrnigsgemaflen Verfahren einfache mafhematische Memodeu angegeben, 
ism MeBwerte entsprechend aufberaten zukdmien unddamit das neuronals Netz trainieren zu kftnnen. 

Besonders vorteilhaft warden Methoden angegeben* urn verrauschte MeBwerte aufeuberdten, bzw; urn 
MeBwerte aufzuberdten, welche eine bekannte and due unbekaunte Rauschkomponence enthalten, urn damit s 
das neuronale Netz auf einfache Weise and effiarienttraimerenzakSnneiL 

Im fofgenden wird die Erfindnng anhand van Figoren wetter erfiutert 

Fig. I zeigt einiekurrentes neuronales Netz 

Fig, 2 gibt ein Beispid des erfindHngsgemfifien Verfahrens an 

Fig.3zeigteineZdtreihe 10 

Fig. 1 grot em Beispid euiesrekurrentenn^ 
mit ein em Systemverhalten SY trainiert, Welches in Fig. 4 nSher dargesteUt ist Am Netz sind freie EmgSnge i 0, 
20 mid freie Aosgfinge 6% 70 vorgesehea Bd Usher angewaudten Trainingsverfahren ffir sofche reknrrenten 

Wdter sind am Netz Rftckfuhiuugen 30 and 40 vorgesehen, wdche bdm rckmrenten Netz ein quasi Gedacht- 
nisreafis£eren,IhFlg.l ist weherbin eine Sdmhtfmie 50 angedeutet, entlang der nach dem erfmdungsgemlBen 
Verfahren die Rfickfflhrungen tew. die Leitungen zwxschen den Ausg&ngen mid EingSngen 30 und 40 aufge- 
trennt werden kdimen. Vorteilhaft wird so dorch das erfmdungsgem&Bc Verfahren der EfFefct erzielt, daB beim 
Training des rekurrenten Netzes audi die Emg&ige 30 imd 40, wdche sonst geschlossen warea, fflr den to 
Thuningsvorgangzur VerfQguag stehen. Dadisrch, daB beim Trainieren des Netzes mehrere Sngfinge verwen- 
det werden, wird ein schneOerer Traiaingserfolg errdcht Bd der Bereitstdinng der TVaimngsdaten fQr diese 
zusitzlkifaenEmgtagewifd^^ 

nacheuumder anfalfenden MeBwerten does technischen Systems beschickt werdm soli tmd dieses nadh dem 
Training nachbfldca sofL Da* so die erfmderische Idee, die beiden zusatadichen EingSnge 30 and 40 ebextfalls 25 
Bestand toil der Zeitreihe von MeBwerten simi, kami fiber die restiichen MeBwerte end deren Fehlerverteflmigs- 
didrte der Wert erredmet werden. Dies geschieht flber das erfindungsgemdSe Verfahren xrrittds der angenom- 
menen Fehlerrcrtefltmgsdichte fiir die EmgSnge 30 und 40 und der Ziemmg von Monte-Carlo-Proben gem£B 

dialer Vefolerm+eatmgsAirht* imrl Mifrtel worrfwlrfimff H» gp^rtgrgnftn Probcn. So ergeben irfdl ZUSafrrifdie TVai- 

nmgseingangswerte gem&B dem erfmdungsgemaBen Verfahren. Die genaueren Details dieser Berechmrogsme* 30 
thode mit dem erfmdungsgemaBen Verfahren sind bd den Erl&utenuigen der Fig. 3 und 4 dargesteUt 

In F|g. 2 wird der erfindungsgemaBe Ansatz nochmals anhand ernes Bdspides verdeutSchl. In Fig. 2 ist 
dasselbereksnrentenraronaieNetz wie inFlg. 1 dargesteUt Bd der erfmdungsgemaBen T itofngsyariante wird 
davon ausgegangen, obwohl das in Realitftt nicht der Fall ist, dafi die beiden dorch Anftremumg der Leitung 50 
entstehenden Eingftnge 30 end 40 mk Werten einer Zehrdhe beschickt werden^ wdebe sfcfa atts den Fingangs- 35 
wertea der Fingange 10 tmd 20 etrecfanen lassetL Per Bngang 30 wird also bdspielsweise mit dem Wert M mid 
d«- FeMerverteSbmg d ausdarF^4bradbdc^w^chermderZeitrememePo 

40 wird dana bdspielswetse mh dem Wert L ans F^g; 4 und dessen Fefaleryerteiiungsdichte g2 brschnckt. Dieser 
nminrtmderZeitr^ediePositm^ 

chen Veizdgenmg der einzelnen Werte innerhalb der Zeitreihe ent^echen. Zl entsporache in dussem Fall einer 40 
Zdtverzagerung vcmyt-i bis yt-2 und Z2 entspridie in diesem Fall an der Ldtung 40 einer Zehtverzdgerung 
zwischen y t ~2 und Mit der erfmdungsgemaBen Methode wird dn Verfahren angegeben, mit weichem 
fehleiKie Werte einer Zeitrdhe mit hcher Genafiigkdt Qber derea Fdrierverteihmgscficfate und 23ehen von 
Monte-Cario-Proben mit anschliefiender Mittelwertmldnng nacbgebildet werden kOmien. Der Gnmdansatz der 
Erfindnng macfat sich diese Idee zunutze, indem die fcflnstKch erzeugten Fingfinge yon rekurrenten neuronaien 49 
Netzen mit kunstlich erzeugten* durch die erfmdimgsgemgfie Methode berecfaneten Wertea, welche aber sehr 
genausicd. besrhtdct werdenu Hierdurch wad vortdmaft eine beasere Tratninssleistmigp welche sich in genaue- 
ren Gewichtenanden einzemen Meuronen wktespiegeh, erztdt Es ist mcht ^forderiich, zur DurchfQhnmg des 
erfinduogsgen^Ben Verfahr ens afle Leitungen, <L k RQckfuhnmgen des neuronaien Netzes, au fai trcnnen, vid- 
raehr kdmien auch faHweise iedigikh einzdne Leitungen auf getrennt werden, da so erheblich weniger Rechen- 50 
anfwand bd der Ziehung der Monte-Cario-I*roben and der Berecfanung des Wertes, flber die Zugrnnddegung 
der Fehlerverteflungsdicfate der restiichen Zehreme, getrieben werden mufiL Bd kompleseren Vorgangen kann 
es auch sinnvoll sein, mehrere IUlckfQhrungsleitungen des neuronaien Nefcces, als hier in Fig. 1 und 2 dargesteflt 
sind, aufzutrennen und auf kompiizierte Wdse die hierfQr erforderfichen Rmgangswerte am bestimmen. VorteH- 
haft wird mit dem erfindungsgemaBen Verfahrai dn Iterationsverfahren angegebent wdches auf Basis der 5$ 
erfmdungsgemaBen. Methode eine Berechmmg auch mehrerer Werte enndgiicht. Faliwdse wird der Fachmann, 
welcher die Erfindnng ausfuhrt, entscheiden kdnnen, welche Methode zur Dundifiihning der Erfindung fQr ihn 
die gOnstigste isL Er wird dabd abwfigen, ob er hinrdchend Rechenkapazitat zur VerfQgung hat, um das Netz 
genau modelfieren zu kdnnen, oder ob es Ami wichtiger ist; einea schnellen Erfolg beon Trainieren des neurona- 
ien Netzes zu errefchen. 60 

Fig. 3 zdgt due Zeitreihe von MeBwerten, welche beispxelsweise einem neuronaien Netz zugef Qhrt werden 
kdnnen. GemfiS ihrer zehlichen Abfolge werden diese MeBwerte beispielswdse von emem technischen System 
erfaBt und gemfifi ihrer zeidichen Abfolge mit yt bis y t -6 bezdclmet Die dargestellten PfeHe zwisclwn den 
einzelnen Kflstchen symbolisieren die Abhiingigkeit^i der verschiedenen Werte tmtereinander. Beispielswdse 
wird in Fig. 3 davon ausgegangen, daB der Wert yt-2 fehlc Die im Markov blanket relevanten Werte, als £5 
benachbarte Werte dieses fehlenden MeBwertes, sind yt-^, yt-s yt-t und yt. Ein solch fehlender MeBwert in 
einer Zeitreihe kann beispidsweise dadurch entstehen, dafi znm fra^ichen Zd^unkt das MeBgerSt, zur Werte- 
aufnahme nicht funkxionierte, oder daB es zwischen emzemen gemessenen Werten gunstig erschemt um das 
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MurcmateNetzbeKCTmtnrim^ 
also noch dem erfndungsgern^n Verfa^ 

usgeffttgen, daB der Wert yt-3 fehlL Die an Markov blanket relevanten Werte, als benachbarteWerte dieses 
fefalenden MeBwertes, sind y t -s, yt-4» yt-2 und yi -4- Beim vcrgesteilten erfadungsgemiBen Verfahren entstc- 
hen diese fehlenden MeBwerte jedoch dun* Anwendimg des crfindungsgcmfiBea Verier* auf retapeme 
neuronafe Netzc» Indqn cfic Rfldkkoppltittgslemingen 30 und 40 a ufg e U emi t werdcn, ergcbcn gen trier pecpieis- 
weise zwei zustaiiche Emg&nge, welche bd Anwendung des erfimiungsgein&Beii Verfahrens zum Turing 
wr^tzl^h nut Lemdateo besdtickt werden konnoL Dies hat nach der erfinderischea Idee zur Fotg^ daB das 
Netz ndt besser zntreffenden GewiAtungsfaktaren ansgestattet wird Db ist der Fall, obwohi dieWote, 
welche flber diese Engfrnge zugeffihrt werdea auf die erfmdungsgemfiBe Art qjwsi ktasdich wenteu 
Dabei ist allgemdn zu bemerken, daB die Anzahl der aufzutrennenden ROckrahrlertungen ledigflcn von der 
Ausprlgung des rckurrentca nearunalen Netzes abh&ogt und vkbt dureh das erfindtmgsgemffle Verfalrren 
beschrfinkt wild. Es kann fellweise andi gfestig sein, racht alle dieser Eingangs-Aiisgai^-Verbirrfuiigwi atifea- 
ttymi^ frrx^ rf« HqIvh m hgtregwnden geringeren ttechenautwandcs eine bimreichende Lemleistungdea 
r^Xetzes erzteDar ist. 

zeigt die Zehrahe aus Ftg>3 in Vcrbmdung mit eraem rekmreoten neuronalen Netz NN, Es ist m 
erkennen, daB y eine zeitabhSngige Variable darstellt, welche das Systemverhaiten SY e rnes tcc hncchen Sy- 
stems rcprfisentierL Wie erkannt weiden kann, entspredien die Werte y» bis yt-6 MeBwerten, wekhe dem 
SystemverfaufSYentnojniiM 
20 wert,dafl diese MeBwerte dem neuror^^ ^ 

Wie audi in Fig. 3 ist der fragliche MeBwert U fflr den Zeitptmkt yt-2 iw*t vorhanden. Fflr diesen MeBwert 
M ist sane Wahrsdiehnichkdtscfidrte el augegeben. Diese Wahrsr hri nl k A k i nt rfidite s kaim begpiebw^ 
nach dem ernndungsgemlBen Verfahren ans einer vorgegebenen bekannten Fehlerverteaungsdicme der ubn- 
gen bekarmten MeBwerte r5ckgcrechnet weiden. Insbesondere wird dabei ansgemitzt; daB sick der fehtende 
25 MeBwert zwischen zwei bekannten MeBwerten befinden mnB 

der benachbarten und der restDchen MeBwerte der Zeitreihe begrenzt wird. Die zngrandetfegende Zeitreihe 
RBtsich wie folgt beschreiberi: 

yt - l&t-uyt-ai.—yt-rf+et (i) 

39 Dabei ist f enlweder bekannt oder wird hinmchend dtirch etn neuronales Netz modelfiert ct bedentet dabei 
eanen additiven unkorreliertea Fehler mit zehiicheni Mittelwert a Dieser Fehler weist dabei — imd das ist ftr 
das erfindungsgemSBe Verfahren essentadl — cine bekamxte odervorgegebcne WahrscheinEcWaitsd^teP«(^ 
auf and versinnbildlidit typischerweise die unmodelHerte DynamSk der Zehreihe. Beispidsweise soil fur eme 

as soldie Zeitreihe, die nach dem oimdiingsgemaBen Vofahren komplettiert werdea soil ein zukunftiger Wert 
vorhergesagt wade*. Dabei bt zu beachtcxv daB zukQnfdge Werte rdativ ai der momentanen gewShlten 
Zeitpontion zu verstehen sind. Das heifit fflr einen Zehpunkt yt -s ist der Zeitptmkt yt-4 can zok flnf tige r Wert 
Unter diesen Vonuissetznngen fiBt stch (fie bedingte WafarsdteinlichkeitadBdilc Mr esnea TOrherzusagenden 
Wert der Zeitreihe wie folgt beschreibea 

Pfcrtjyt-uyt-zi-^yt-ii) - PeCy- Hy*~ uyt~*-* n-ni) (?) 

Wie bereits erwfthnt mnB die FdilerverteUungsdidite bekaimt sent Diese Verteimngsdidite k ann en tweder 
anhand des Systemverhaltens und bekaimter anderer &nfierar GrOBen ennit^t oder vorgegeben werdea Bine 
typische Fehlerverteflung. die in der Praxis auftiitt ist de GauBverteihing. MH einer solchen angenommenen 
GauB r sdieii Fehlerverteihing IfiBt ach die bedingte Wahrsdieinlichkeitsdidite wie folgt beschreibem 

PCyt|yt-i»yk-2.. , yt-w) -G^t^-i—^yt-Nio*) ($> 

Darin bedeutetG^c^o 2 ) die Notation fBr eme nonmleDicite, die beixbestimmt wird n^ 
und einer Varianz o 2 . Geht man davon aus, daB das zu bescfareibende System in Form einer Fol^ von Werten 
auf emer Zeitachse dargestellt wird* so kann man die einzemea W«te von yt auch ats Zufallsvariable in euiem 
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he vorherzosagen, indem die vorhandene Information aus den restlichen Werten mdgCchst voDstfindig verwen- 
det wird Unter Voraussetzung der Annahmen, die zuvor gemadit wurden, l&Bt sidi die gesamte Watascheinfich- 
keitsdichte der Zeitreihe wie f cigt beschreiben: 
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Dabei wird davon ausgegangen, daB yt^k mit k £ Nderfehlende Wert isL Voransgesetztp daB y* = {y t -k}und 
_ \y t _ u „, yt-k-NMyt-rf gih, kann der erwartete Wert der in der Zeitreihe vorherzusagen ist wie folgt 
beschriebenwerden: 

ECyt|Mt-i) - Jlt>t-!,-.-,yt-^^y«-N)PtyV^ (5) 

Dabei gehen folgende Voraussetzungen: 
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Mt-i stefat f&r afle Messungea bis zum Zeftpnakt t— L Die voranstehende Gieichung ist die gnmdlegende 
Oeichong {Or die Vorhersage nit fehlenden Daten. Dabei ist besonders zu beachten, da& die Unbekannle yt-k 
nicht nur von den Werten der Zeitreihe vor dem Zeitpunkt t — kabhfingt sondera audi von. den Messungea 
irach t — k. Der Gxund bestebt darin, daft die Variablen in y™Uyt ein mintmqlpc Markov blanket von y*_k 
fonnen. Dieses minimaie Markov blanks bestebt was den direkten Vorfabren and den direkten Nachfahren 
einer Variable und alien direkten Vorfabren von Variable des direkten Nachfolgers. Im betracfateten Beispiel in 
Fig, 4 aind die direkten Nachfahren y t . . . yt -k + 1- Di e direkten Vorfabren sind: 

yt-k-i...yc-k-N 

nnd die direklen Eltem des Nachfolgers der Variablen aind: 
yt-i.-*yt-k-M+i- 

Ans den theorerischcn Grandkgen ist bekannt daB eine Variable imah hm gig von einer anderen Variablen is 
dieses Netzwerkes ist wenn die Variablen mnerhalb des Markov blankets bekannt and Deshalb wird die 
ben6tigte becfingte Dtdite ans Gleidning (5) wie folgt bestimmt: 

PCs^y"*) t» P6rt-ifartwa-.,yt-fe..,yk-i-w) x F(yt-2iyt-3, . ♦ , yt-t - . ^ yt-2-N) . . - P(yt-k|yt-k- i, . y t — 
k-NX (5b) » 

Der hier bescbriebene Pall ernes fehlecden MeBwertes kann aucfa ohne Beschrankung der Erfindung auf 
mehrere nebeneinanrier liegende fehlende Mefiwerte ausgedebnt werden. Fails dies da- Fall ist, mnfi beispiels- 
weise nacb dem erfindiiixgsgem&fien Verfahrea zunSchst der eine Wert anband seiner Nachbarn and Eltern and 
Vorfabren bestiixunt werden imd ink diesem bestimmten Wert dann der wcitere Wot ermittelt werden. Dies 25 
kann solange hin- and hergeben bis eine hinradtende Genauigkeit erreicbt werden. FOr diesen Fall glf= 

C {yt-t,yt~a--*yt-4 (5c) 

FQr afle fehlenden Werte der Zeitreihe zwischen dem Zeitpunkt t— 1 and t — N, wobei wdterhm gilt: 30 
^efyt-hye-a.^yi} (5dX 

welche die Zahl aller MeBwerte fusznmZeitpimktt— 1 reprdsentiert Audi gilt 

Wy")- P&t-i^-.yt.yi) (5eX 35 

wobei <Se recbte Sdte in (5e) ans Gieichung (4) erhalten wird. im aUgemcinen sind diese Integraie in dm 
voranstebenden Gleidmngen fOr die Funkfrra IQ> foils <£es eine nicbtBneare Fonktion ist nicht analytisch 16sbar. 
Details ffir die nmnerische LOsimg mit HDUfe der Ziehmig von Monte-Carfo-Proben werden im Zosammenhaiig 40 
mit Fig, 3 angegeben. FOr den Fait daB ein wetterer MeBwert der Zeitreihe nachgebildet werden soil, wie dies 
bei zwei aufgctrennten Rflckkopplongsvci bindimgen der Fall ist; sieht das erfmdungsgein§Be Verfahren eine 
iterative Approximation der fehlenden Werte vor. Beispielsweise ist far das Training des Netzes zusStzfich der 
Wert L fflr den Zeitptuikt y t - 3 nachzubilden. Ffir diesen MeBwert M 1st seine WahrschehtHchkeitsdichte e2 
angegeben. Diese Wahrschdniidtkeitsdichte c2 kamt beispiebwcise nach dem erfindongsgemftfien Verfahren 45 
aus einer vorgegebenen bekannten Fehlerverteilungsdichte der flbrigen bekannten Mefiwerte rttckgeredmet 
werden. FQr die Berechnung von zwei sofchen Werten L imd M wird zunSchst L beispielsweise als bekannt 
vorausgeseizt oder geschatzt Darans wird M beredmet imd mit dem bekannten Wert M wird anschKeBend 
erneut L bestimmt. Dieser Iterationsvorgang idnft vorzugsweise so lange ab, bis eiw* hmi whpnrfA fw^namgirgHt 
der Werte gegeben ist oder bis das Netzgenaugenugtrainiertist. 50 

FOr den Fall, daB y u yt magliche Werte der Zeitreihe darken, sollra y™ Q Jyi,.. v y t }afle MeBwerte 
bezeicbnen undy^ — {yu * * * yjy 1 ^ aJle unbekannten Werte bezeichnen. Das neuroiiale Netz NN, welches {fie 
Funktion f modeOiert werde betspidsweise mtt einom Satz von Gewichlen w parametrisiert Dann^It: 

f(yt-i.-^yt-N)» NNw(yt-b.^yt-N) ss 

Die logarithmische WahncheinlichkeHsfunktion lautet damu 0 

L - logJRMCy^yt-i,. ^y^ijdy 11 

wobei dann diegoneinsame Wahrscheinlichkeitsdichte sichzn 
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approximiert und f Or das neuronale Netz folgender Zusammenhang (Or die Berechnung der FehlerverteBungs- 
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dichtegQt: 

P"(yt}yt-i,yt-a--,yt-N) -PJfa- NN^yt-uyt-».-yt-w) C 7 ) 

Rlr das I^menmitHafe\^Badc}m)pagatk»n, oder aniierer 
noch der Gradient dcr logarithmischcn WahrschemlicWteitsftinktion benotigt. welchcr ach zuz 

er «ibt.EsHtanauneikeii f daBliierDd 

den FaS, daB one Ganfivcrteitang ftr die Fdblerverteihingw^g^ai^^darmw 

£ f (y, -NN.( jH „w) 3NN ^ (y -- :? "' ) p"y<° . yw* <ea> 

Hjiftfaflg nllr.yr y -Mfaekmirtsm&das Integral vcrofawindet ... * , . 

Falls die McBwertc von emcm jpsltTtidiffn aber bekanntea Rauscfaen Hberiagcrt werden ergcfaea ach die 
fblgenden Zusammeuhfinge. Beispidsweise gilt wieden 

yi - ttyt-i,yt-i- .*yt-M> + d 

In dieser Variante der Erfindung soil jcdoch kcrfn dirdrtcr Zugriff auf yt besteben. Anstattdessen wird die 
Zeitreihe 

«t-yt + 8i 

ranraxcn. Darin brdnrtrt St ^hhi^^ JtansAen nut Mittclwert NiilL Unter der Voraussctzung, daB z ^ 
{zi^.zt-tftnriy = {yi~.ytjgelten,^ 

fc=ttH. W 

DamUlaBt rich die Redxawor«&riftfflr^ 

E(yt|z - ffi(yt-i,-. ,5t-N)P(yt-i.-.,yt-Nlz)dyt-i ...dy t -w (9) 

ir. m Hbt Otfiem der Wahi^Hnnrhkrifsftrnktimi ffir das Training benedmet widen. Pgr den Fall, 
daB eme GauBverteihmg des Ranschens nut 

vx*riiegt,er-gibtsich: 

Z J(y,-NN w Cy M .„,y t .H)) J£ — 

Dem nenronaien Nets werden ba einer Variante des erfmdangsgcmiBen Verfuirens , iHa^ndawefae Werte 
zageffihrt, die verrauscht oder mcbt genan bestimmbar and. Dtnch die Approximation der - Gewwhte aa 
neuronalen Netz werden dabei liber die Fanktkm f. wekhe dibei durcfa das neuronale Nete nachgebildet wird. 
neue Werte der Zeitreihe bestimmbar. Diese neuen Werte der Zeitreihe werden nn AnscbluB dem neuronateii 
Nets NN zngefQbrt, welches damns wiederum durch Nacfabudung der Ftmktion fneue Werte der Zertre&e 
bestnnmt. Dieser iterative Vorgang wird solange fortgesetzt, bis eine hinreicbende Groaragkeit der zu bestnn- 

"aJ^^Be^w^Wilender Werte nut Hilfe der Monte-Carlo-Methode wird von folgenden Grundla- 
gen ausgegangen. Es ist hier za beachten, daB aHe Ldsnngen die Form 

p<n,m)P(u|in)du (9b) 
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aufweisen, wobei u den Satz von unbekannten Vanablen und m den Satz von bekamrten Vanablen bedeutet Ein 
Integral dieser Form kana beispielsweise gelfist werden, mdem Zandfaproben der unbekannten Vanablen 
gem&B P(o|m) gezogen werden- Beispielsweise warden diese Proben mh a 1 ,., u* bezetcimet Darans ergibt sich 

fnlgt*ptfer y^yfttfiin<*nfaaitglftrdfe Afmflfienmgi 
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Es istzu beachten, dafl m dieser Glddiung u den Wert y t -k, wdcher feblt, entspricht Mh dleser erfmdnngsge- 
m£8en Losung reduziert sich das Problem also darauf, au$ P(u|m) Proben zu Ziehen. Ftir den Fall, daB ledEgficb 
eine Variable fehlt, also beispielsweise ledigEch eine Rfirkfflhnmg aufgetrennt wurde, reduziert sich das Problem 
also auf das Probenziehen aus etner einvariablen Verteilung, wetche mit Hflfe des ^sampfing4iiqxirtance-r«am- 
pfing^ oder anderen sampfing-Techiiikeii [1] getan werden kann. 15 

FQr den Fall daB mehr als em MeSwert fehlt, wird die Situation etwas fcomplfederten. Der Grand besteht darin, 
daB die unbekannten Vanablen in der Regel voneinander abhflngen und daB von der Vertdhmg aller upbekann- 
ten Vanablen gezogen werden moJL Eine aflgemeine Losung dafQr g&t das Gibbs-Sampling [t] an, Beim 
Gibbs-Sampling werden die unbe kannte n Vanablen mh Zufaflswerten, oder besser mh geschiftzten Ansgangs- 
werten inidafisiert, weldie beispielsweise aus den benachbarten Wcrten der fehlenden Werte abgelehet werden 20 
fcdnneo. Jm Anschlnfi wird due der unbekannten Vanablen tq ausgewtfah and eine Probe von P(ui[nyi/uj) 
gezogen; dannwird us auf diesen Wert gesetztDaraufhin wird die Prozedurfttr die nSchstc unbekannte Variable 
wtederholt osw. Abgesehen von beispielsweise den ersten Proben werden die Proben vorzugsweise mit der 
korrekten Febkrverteifan gsdfchte gezogen. Dies bedeutet jedoch, daB fllr afle Unbekamrtgn, ale jemafa m der 
Zeitreihe auftraten, Proben gezogen wertlen mflssen. In der Praxis kann aber beispielsweise das Zehfenster aus 2s 
detn Proben gezogen werden auf eine vernfinftige GrdSe beschrinkt werden. Beispielsweise karni cfiese Grdfie 
der GrSBe des Markov blankets fDr die fehlenden Werte entsprechen. Hierbei ist zu beachten, daB fflr den Fall* 
daB zwischen zwd fehlenden Werten N anfemanderfolgende Werte bekarmt sind, die Kopplung zwischen den 
Unbekannten aufbricht end daB weitere Werte der Zeitreihe deshalb nicht berffckrachrigt werden nxQssen. 

Das Probenziehen fur zuicGnftige Werte ist nach dem er6ndungsgemaBen Ver&hren besonders einfach. JBs ist so 
jedoch za beachten, daB es so rriefct fOr detenrnnisrisrhc Systeme funkdoniert Durcfa die erfmdnngsgemABe 
Vorgehensweise wird eine besonders einfache L5scng fOr diese besonders kompfiziert erscheinenden Sachver- 
halte gefunden. Beim Training des Netzes wird beispielsweise der Durcbschnitt der FeUergradienten gebildet, 
mdem zu deren Berechnnng die Werte der Zeitreihe verwendet werden, wekhe mit den Proben bestimmt 
warden. Beim Ziehen der Proben nach der Monte-Cario-Methode kann beispielsweise wie folgt vorgegangen ss 
werden. 

Es sollen beispielsweise K-$chritte in der ^doinft der Zeitreihe vorhergesagt werden. Im Zusammenhang mh 
don oben besproclien» Gfeichnngen bedeutet dies„ daB die Werte yt—u * yt-R+i fehlen, und daB imtc 

di«Ren VnraMmftMtwgftn wrh^rgwTOgf wectten mil. Ilirtar rfiAgtm Vrmn«<^imgi«ii nt Ai+ MoCte-Carto-Me- 

thode sehr emfadL Zuntcbst nmB beispidswdse von der VerteSung F(yt-tc+ijyt-K. . , yt-K-N) eine Probe 40 
gezogen werden. Diese Probe wird mh y*t-K+i bezeidinitf, Mh dieser Probe and den vorangegangenen 
Messungen wird beispielswe^e eine weitere Probe yt-K+2 aus der VmeOung P(yt-K«|y , t-K+ L-^yt^K+i-N 
usw. bis jede Probe fOr jede Unbekannte erzeugt wurde. Wenn man diese Prozedur S-mal wiederbolt, so erh&it 
man 



Die experimenteOen Befunde haben gezeigt, daB mit dem erfindungsgeica&en Verfahren gute Trainingsld- 
stungen erzieit werden. 

Lheratur 



45 



so 
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Patentansprdche 

1. Leiuveifahren furdn reknrrentes neuronales Netz mitfolgenden Merkroalen: 

a) mindestens erne von ehiem Ansgang za dnem Fjn g gn g zurQckgenlhrte Verbindungsieitung des 
N^zes wird zum Lemea aufgetrennt and der entstehende ^" g«"g in Form eines Lerneingangs 



es 
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b) dieEing^ 

fa fit* 

c) von dnem den EmgangsgrfiBen Qberiagerten unkorrefierten Rauschen endGcter Variant das bn 

d) dkd«n Umeingang zugefDhrte HmgangsgrtJBc wird *k feblender MeBwert " d« Zeilrefo m 
Form ernes Fehlwertes aufberatet, mdem mmdestens aos dem Fehiwert m dwZwtfeaebcnacMbarten 
MeBwertenund der bekannten staQstischen Raasdwerteifang dessen statistrache Jduwert-RauscJwer- 
tefflang berecfanet wild end fur den Fehiwert gemaB der FeUgert-Ranadwerteatang m m dn h aw erne 
Monte-Carfo-Probe des FebJwertes gezogen wird, die an seine Stelle in der Zota^e tntt 

2. Verfahren nach Ansproch 1, bei dem ein vorhenroagender Wert m der Zeitrribe ab i ™S n ™wt 
bestimmt wird, mdem fOr den Fehiwert menrere Monte-Carky-Proben g^^^ t ^ l ^nosewerte 
bestimimwi^dWvrobdalaProgmMerwertdasariA 



a vSatolm^ Ansproch 1 oder 2; bei dem mmdestens ^T^^t^ 8 ^^^^^ 
zarOdcgefOhrte Vcrbindungriatungen des Netzes mm Lernen anfeetremrt .werden 
tfagSpml^emaefStoa^ 

arfdieseWeiseenraigtenzwei fcMendenoiiduimuttcIbar bei 

SSd ^^Seingang born Training zugtfuhrt werden, znnachst der Erste arfberertet wird and 
im AnscntaB dar*n nnter Zohflfenahmc des zuent Mifbereheten der Zw«te bestimmt wird. 

4. Verfahrea nach einem dervorai»gebeadenAnspr&Ae,WP»rhrt m e hrfarfi ansgrfOhrtwn-d. 

5. Verfahren nach emem der vo«iiiEelien^ 

garion die l^njsdmttweite fflr die auf 1 normierten EingangsgrOBendes nearonalen Netzes zu 04 drwdiert 

GauBverteihmgfestgdegtw^ , . , _ 

7. Veifahrcn nach einem dervorangeftenden Ansprflche nut oner Zertraihe der Form: 

yt - *t>t~uyt-*"*yt-N)+et (1) 

niitietistatistische Rauscliverte3iing 
y: MeBwerte der Zeitreihe 
y t : vom neuronalen Netz TOrfoerzusagender Wert 



Vt:vom neuronalen Nctz TOrfoerziisagOTder Wert b ~ r . 

Wd«dieI*iifclionf 
tezu: 

P«(yk-i - *Xy*-i,yt-*- ^*-H» - P(ytlyt-uyt-2-.,yt-N> (2) 

bestimmt wird. woraus die Moiite-Cario-Probeii y,-* 1 , - • » yt-k* gezogen werden tind aich ein erwarteter 
vom neuronalen Netz vorfaerzasigender Wert mit den gezogenen Proben zu: 

ergibt b 
mit:yt-~k feblender MeBwert in der Zeitreihe 

mi alte bekannten MeBwerte der Zeitreihe 

?VeS^ mcbSem der vorangehenden Ansprifcte beidemdasNeuronaie Netz mit mindesteiis emem 
arfbereitete n Wert nach folgender Lernfunktiott trainiert wird: 

BL 

mit: w: Neuixmengewicht 

L: logariinmisdhe Wahracbeinlirfiketofunktion 

n: Lernfaktor 

wobeigilt: 
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mit: NN W : Werte der Funktion aus dem neuronalettNetz 

wcfaei fdr y*i McBwerte der Zeitreihe Verwendung fmden und, foils ein Wert w±t vorhanden ist; aus der 
Wahrscfaeinfichkeitsveiteilangsdicfate: 

p»w> s 

Monte-Cario- Proben gezogen werden mit: 
mi alle bekanoten MeBwerte der Zeitreihe. 

9. Verfahren nach Anspnich 6, beidem das Neuronale Netz mit nrindestens einem aafberatetea Wertnach 
fblgender Lerofunktkm trainiert wnrcfc io 

dL 

15 

mit: w: Neuronengewicfct 

L: Zogarithmische Wahi^einlichkehsfunktioa 

Tj:Lernfaktar 

wobeigih: 

20 

dw S^jSi 

25 

nut: NN«: Werte der Ftmfction aus dean neuronalen Netz 

wobei ffir y*i Meflwerte der Zeitreihe Verweodtwg findea and, falls ein Wot nfcht vorhanden ist, aus der 
Wahreeheinlichkgitsverteg 

Monte-Garto-Proben gezogeo werden mit: 
m: alle bekannten MeBwerte der Zeitreihe. 
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